Functional approach to the 
fermionic Casimir effect 



C. D. Fosco and E. L. Losada 

Centro Atomico Bariloche and Instituto Balseiro 
Comisidn Nacional de Energia Atomica 
R8402AGP S. C. de Bariloche, Argentina. 

Abstract 

We use a functional approach to calculate the Casimir energy due 
to Dirac fields in interaction with thin, flat, parallel walls, which im- 
plement imperfect bag-like boundary conditions. These are simulated 
by the introduction of 5-like interactions with the walls. We show 
that, with a proper choice for the corresponding coupling constants, 
bag-model boundary condition are properly implemented. 

We obtain explicit expressions for the energies in 1 -|- 1 and 3-1-1 
dimensions, for massless and massive fields. 

1 Introduction 

In recent years, a renewed interest in the Casimir effect [1] emerged, as a 
consequence of new, refined experimental techniques [2], which raised the 
required standards of the theoretical calculations. They were usually based 
on gross simplifying assumptions, mostly about the properties of the mirrors. 
To explain the experimental results, it is then important to use more accurate 
models, including corrections due, for example, to the imperfect nature of the 
mirrors. 

In this paper we are concerned with the calculation of the fermionic 
Casimir effect for imperfect 'mirrors', in a functional integral framework. In 
a fermionic context, the mirrors are meant to partially refiect the fermionic 
current; including as a particular case the 'bag' boundary conditions, where 
the current is completely reflected. 

From a technical point of view, we shall use an adapted version of a 
previously used functional approach [3], to cope with the case of fermionic 
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fields and imperfect boundaries. In the fermionic context, this method had 
already been applied to the case of fermions in 2 + 1 dimensions satisfying 
bag conditions on a curve and coupled to a gauge field [4]. 

To cope with imperfect mirrors, which will occupy parallel planes, we 
follow [5] to introduce a coupling of the bilinear '^{J^|: to a 5-function potential 
of the proper coordinate. We shall then calculate the resulting Casimir energy 
as a function of the distance between two mirrors, using to that effect a 
functional formalism which maps the original problem to a one-dimensional 
one, where the fields live on the surfaces of the mirrors. We show that 
exact results may be obtained by for 1 + 1 and 3 + 1 dimensions; for the 
1 + 1 dimensional case our result in the perfect bag case agrees with the 
one of [5], in spite of the fact that the treatment of the problem is quite 
different; in particular, the description of the interaction cannot be mapped 
in a straightforward manner from one calculation to the other. 

This article is organized as follows: in section [2] we introduce the method 
in general; in [3] we apply it to the calculation of the Casimir energy for two 
imperfect mirrors in 1 + 1 and 3 + 1 spacetime dimensions. In section H] we 
analyze the relationship between this kind of interaction and the bag bound- 
ary conditions, which emerge when the coupling constant takes a particular 
value. Finally, in section [5] we present our conclusions. 



2 The method 

In this section, we shall introduce the method used for the calculation of 
the vacuum energy, in the presence of 'defects' which generate the approxi- 
mate bag-like boundary conditions, for the case of a fermionic field in -D + 1 
spacetime dimensions. 

Our approach is based on the well-known property that the vacuum en- 
ergy, Eq, may be obtained from the Euclidean functional Z = , as: 

r 

^0 = lim - , (1) 

where T is the extension of the (imaginary) time coordinate. On the other 
hand, in the functional integral representation, Z is given, for a Dirac field, 
by: 

Z = [ Vi)V^e-^, (2) 
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where S, the Euclidean action, will be such that S = Sq + Sj, with 



Si 



'0 



/ 
/ 



X iIj{x) + m) ipi^x) 



X ipi^x) V{x) i/j{x) 



(3) 



which are the free and interaction actions, respectively. 

In the last equation, V{x) denotes a 'potential', introduced in order to 
simulate the boundary conditions. Although we shall mostly consider just the 
cases of one or two defects, we write (for the sake of generality) its expression 
for a system with equally-spaced defects: 



where L and a = 0,l,...,A^ — 1. ^fis the coupling constant, whose 

role in imposing an approximate boundary condition is discussed in section [H 
It is worth noting here that, as mentioned in [5], Dirac's equation with this 
kind of potential is not well-defined. However, also in [5], a consistent cal- 
culation of the Casimir energy for a Dirac field in 1 + 1 dimensions could 
still be performed, since one can set up the problem in terms that avoid the 
calculation of the eigenstates of the Dirac Hamiltonian. We shall here use 
an approach that also bypasses the eigenproblem for the Dirac Hamiltonian, 
relying instead on the corresponding Euclidean propagator. The reason this 
can be done is, of course, that the vacuum energy may be expressed as a 
function of the expectation values of fermionic bilinears which are, at least 
in this case, well-defined even when the eigensystem is not. 

Of course, the total energy for this kind of potential will be, since there 
is translation invariance along the coordinates parallel to the surface of the 
mirrors, proportional to the 'area' of each defect; thus we will be interested in 
evaluating energy densities, or pressures, rather than the extensive quantity 
Eq (except the D = 1 case where the two quantities coincide). 

In D spatial dimensions, the defects are co-dimension 1 hyperplanes. 
Putting the system in a d = D — 1-dimensional 'box' of side a, we obtain the 
vacuum energy density So dividing Eq (calculated for such a box) by a*^, and 
then taking the a — >■ oo limit: 



We have denoted by x the coordinates in D + 1 dimensional spacetime x = 
(xo, . . . ,xd)- It will turn out to be convenient to introduce also coordinates 



N-l 




(4) 
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y G M'^"''^, for the subset of the d+1 spacetime components that parametrize 
the a*"^ defect's world- volume, as a; = {y,aa)- Analogously, momentum- 
space coordinates in D + 1 dimensions are denoted by p, while q is reserved 
for its d + 1-dimensional restriction: 

P = {q,Pd) , (6) 

where is the momentum component along the normal direction to the 
mirrors. 

In terms of these conventions, we have: 

N-l 

Si _ TT - g f d'^+'^y 4){y,aa)^{y,aa) 



e 

a=0 



JJ^ g-9 '/'(y.^c) V'(2/.aa) (j-j 



which, following the methods of [3], may be written in an equivalent way by 
introducing a pair of auxiliary Grassmann fields Xa{y) and Xaiy) for each 
defect. The vacuum functional then adopts the form: 



Z = Zo [ [YlVxaVXa] exp[j2 f d'+'yxl{y) 

a=0 a=0 

+ g f d^+^x / d'^+^x' f]{x)G^''\x,x')r]{x') \ , (8) 



where: 

while Zq is the vacuum functional in the absence of defects: 



(10) 



Zo = j ViIjViI) e^^° = det ( ^ + m) . 
On the other hand, 

G^'^^x^x') = {x\{^ + m)-^\x') , (11) 

is the free fermion propagator. Note that the auxiliary fields are naturally 
defined on each one of the mirrors' world- volumes. 

Since we are interested in calculating the dependence of Sq on the posi- 
tions, a^, of the defects, we shall discard any constant which is independent 
of those variables. Constant factors in the functional integral, like Zq, gener- 
ate precisely that kind of constant, which we shall therefore ignore (without 
bothering to rename the corresponding subtracted quantity). 
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Replacing (l9|) into ([8]), we may write the latter in the more explicit form: 




7V-1 



Z 



(12) 



where we introduced the matrix kernel: 



(13) 



where each a, (3 element is also a matrix in Dirac space. 
The formal result of evaluating the integral above is then: 



Z = det lCaf}{y,y') 



(14) 



where the determinant is meant to affect all, continuum and discrete, vari- 
ables. This means that: 



Since there is invariance under translations in the y coordinates, we may 
Fourier transform with respect to them. The kernel becomes then block 
diagonal in the continuous variables; besides, a proper counting of the modes 
in the box shows that the a'^T factor is cancelled by an identical one coming 
from the numerator; thus, in the limit. 



where the tilde has been used to denote Fourier transformation, and 'tr' 
denotes trace over the a, (3 and Dirac indices only. 

In the following section, we deal with the evaluation of the previous ex- 
pression for £q in the most relevant case, i.e., N = 2. 

3 Casimir energy 

Besides the subtraction of the vacuum energy in the absence of the defects, 
there is another constant to get rid of: the self-energy of the defects. It may 
be identified as the result of evaluating for ^ — ^ oo. Thus the prescrition 

^Again, a finite spacetime box is introduced before tacking the limit. 




(15) 



where 'Tr' is the trace over continuous and discrete indices Q. 




(16) 
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for such a subtraction, which will render a finite quantity as a result, is to 
subtract from 80 its limit when L ^ 00. Of course, to give meaning to the 
expression at intermediate steps, a regularization is introduced. 

Since the last subtraction brings into play the self-energy of the defects, we 
present first = 1, where the self-energy first emerges, and then the N = 2 
case where we evaluate the energy for the case of two defects, identifying and 
subtracting the contributions of the self-energies of the two mirrors. 

3.1 One defect 

For = 1, the K. matrix has only one element (a = /3 = 0), thus there is no 
need to use a, (3 indices: 

m + i ^ 



)C = 1 
so that So becomes: 

£0 = 



9 



Wiq) , W{q) ^ r^-^ 



-trjln [1 + ^Wiq)]} . 



(17) 



(18) 



When D = 1, the Euclidean gamma matrices have the single eigenvalues 
±1. Then the vacuum persistence amplitude is: 



So = - /^Infl + ^ + ^i^ 



(19) 



which has a logarithmic UV divergence, easily regularized by using a fre- 
quency cutoff. 

A lengthier, but quite straightforward calculation shows that, in D = 3: 



i2n) 



In 



4 



gm 



(20) 



which is quadratically divergent. Again, an Euclidean can be used to give 
meaning to the integral, whose explicit form we do not need: indeed, it will 
emerge in the next subsection only to be subtracted in order to fix the energy 
to zero when L ^ 00. 



3.2 Two defects 

When N = 2 the matrix ]C{q) is given by 

1 + IW{q) f (W{q) - ^d) e-^ 



1 + lW{q) 



(21) 
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with W as defined in (fTTj) . 
Using the property: 



det }C{q) = det /Cqo det /Cn det 



'10 



(22) 



for the determinant of the matrix in (|2T|1 . we may identify the first two 
factors above as yielding the self-energies for the defects, as studied in I3.1[ 
We first introduce an UV cutoff to give meaning to those two factors (the 
third factor is convergent) before subtracting the corresponding self-energies. 
The resulting pressure is: 



(27r)' 



jtT[\nM{q)] 



(23) 



with 



M{q) =1- /CoV/Coi/Cr//Cio . (24) 

Let us now evaluate Sq above for D = 1 and D = 3. In one spatial 
dimension, if we take into account Dirac's indices, /Cq,/? is a matrix with two 
2x2 blocks, while M is just a 2 x 2 matrix. A somewhat lengthy but 
otherwise straightforward calculation shows that the latter has eigenvalues 
Ai, A2 given by: 



Ai 
A2 



1 

1 + 



16 ql 



-2L|<7o| 



(25) 



[(4 + g'^)u{qo) + Agmf 

where uJ^qo) = q^ + w?. Therefore, the trace in the expression for the 
vacuum energy can be calculated exactly. In the massless case, we can obtain 
the exact expression for the energy as a function of g^ since: 



TT 



dqo In 



16 g' 



-2Lqo 



(4 + ^ 



2^2 



or, 



£0 



where 



-—Li J- 



-16g' 



(26) 



(27) 



(28) 



2 7r~"^L(4 + ^2)2J ' 

is a dimensionless constant which determines the strength of the Casimir 
interaction, and Li„(x) denotes the polylogarithm function. In figure ([I]), 
we plot the dimensionless combination Ci{g) = —L£q as a function of the 
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(also dimensionless) variable g, for the massless case, and for different values 
of M = rriL, corresponding to the massive case (where there is no closed 
expression for the energy). Note that there is a maximum at g = 2. That 
value, as explained in the next section, corresponds to exact bag boundary 
conditions, where the boundaries are more effective and as a consequence the 
energy reaches it maximum possible value: 



IT 



g=2 



24L 



(29) 



Note that in the approach of reference |5j, bag conditions correspond to an 
infinite coupling constant. The difference is explained, of course, by the 
different meaning of the coupling constants in both cases. For example, the 
propagator in our approach includes tadpole-like propagation between the 
wall and itself, which naturally change the effective value of the coupling 
constant. Nevertheless, our g = 2 agrees with A ^ oo. 

A fuller explanation of this fact is given in section [H where we derive the 
form of the fermion propagator in the presence of the defects, studying its 
behaviour close to the boundary. As already mentioned, for the massive case 



0.12 



M = 




Figure 1: Ci{g) as a function of g, for different values of M = mL. 



the corresponding integral cannot be performed exactly. In Figure 2 we plot 
the result of numerically evaluating that integral for Ci{g,mL) = \Sq\ L as a 
function of M = mL, for different values of^Q'. 

In 3 spatial dimensions, each block in is a 4 x 4 matrix. For a massive 



8 




0.0 0.2 0.4 0.6 0.8 



Figure 2: Ci as a function of M = rriL] D = 1. 



field, and using the notation u{q) = i/m^T^ find that: 

X:ooXi/C^//Cio = — — — ^ -2 -2q^X + gq^-iD 

[(4 + + 4^m] L 

+ i{2m + guj{q)) 4 + i{2uj{q) + gm) ,(30) 

with g = (po,Pi,P2), Pd = Ps and 7^ = 73. 

The exact eigenvalues of Ai can be found after some algebra, the result 
being: 



Ai 
A2 



2 ^2 



16 g^ q 



[(4 + g^)uj{q) + 4gm\ 



-2oj{q)L 



(31) 



each one with multiplicity equal to two. Then, after integrating out the 
angular variables, the Casimir energy per unit area may be written as follows: 



where 

C3(^,M) 



£o{g,m,L) 



dr In < 1 



C3{g,mL) 
L3 



[AgM + (4 + ^2) + NP\ 
9 



2 ( ' 



(32) 



(33) 



which converges in both the IR and UV regions. 

It is convenient to introduce first the result corresponding to the case 
m = and g = 2, where the energy reaches a maximum: 

77r2 



^o(2,0,L) = 

in agreement with [6] (see also [7]). 
When m 7^ 0, and introducing M 



dr r In 



2880 L3 



mL ^ 0, 
16 g'^ e' 



(34) 



(35) 



[igM + (4 + g'^)Vr^ + ]2J 

This integral cannot be found analytically, but it has a convenient form for 
performing it numerically. In [3l we plot c^{g,M) = \Sq\ as a function of 
g for different values of M = mL between 0.1 and 1. We see the expected 
decreasing behaviour of the energy when m grows, regardless of the value of 
g. This happens because the fermionic field propagator has a faster decay 
with distance the larger is the mass, reducing the interaction energy between 
them. On the other hand, the energy always has a maximum when g = 2, 
for any value of the mass. As we mentioned for the D = 1 case, this comes 
from the bag boundary conditions, which are met precisely at this value. In 



M = o 



C3 



0.020 



0.015 



0.010 



0.005 




Figure 3: C3 as a function of g^ for different values of M. 



figure m we present the complementary view, by plotting C3 as a function of 
M, for different values of g between 0.1 to 3. It can be seen that the energy 
becomes negligible when M ~ 2 which means that the distance between the 
plates is of the order of the decay length of the propagator. 
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0.010 



0.5 1.0 1.5 M 

Figure 4: C3 as a function of M = mL, for different values of g. 



4 The (5- interact ion and bag boundary condi- 
tions 

Bag boundary conditions are usually formulated in terms of the ferniionic 
propagator G{x,x') = {ip{x)'ip{x')) , in the presence of boundaries. For the 
case of just one defect located at Xd = 0, one says that it imposes bag 
conditions on its two, 'right' and 'left' faces, when: 



lim (1 ± 7z))G'(x,x') = 



(36) 



respectively. They imply that the vacuum expectation value of the normal 
component of the fermionic current vanishes when approaching the wall from 
each side. Similar conditions may be imposed, of course, on more than one 
wall. 

Following an analogous procedure to the one use for the calculation of Z, 
but now for the propagator: 



{ilj{x)ij{x')) 



J Vi)Vi)ip{x)tlj{x') 



we obtain: 



G{x,x') = G^°\x,x') + T{x,x') , 



(37) 
(38) 
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where: 



T{x,x') = -gY^j ^^e-'^^y-y'^&'\q,xn,a^) 

X [k~\p{q)G(^){q,ap,x'j,). (39) 

with 

G^^\q,XD.x'o) = /^e^^^'^("°-"i')G'(°)(a;,a;') • (40) 
J 27r 

The matrix elements of the inverse of K may be written in terms of the 
matrix elements of K: 

(Cq^ -[ACoo]~^A^oiC'r^ \ 

1 ' ^^^^ 

where [/Coo]~^ and [/Cn]^^ denote the inverses of the respective matrix ele- 
ments (not to be confused with the matrix elements of the inverse) . We have 
used the definitions Co = /Coo — Ai;oi[Ai^ii]~^Ai!io and C\ = /Cn — /Clo[/Coo]~^/Col• 
Let us apply the formulae above, for the sake of simplicity, to the case of 
massless fermions and only one defect in D = 1, we see that, when xi — > 0+ 
and x[ > we obtain: 

(l±7i)[G(g,0,x;) + G{q,0,x[)] = l'"^'^! (4^ ^ ^) 



X 



J =F isg{q)jo ± 7i + «sg(g)7o7i , (42) 



where sg is the sign function. Since we are interested in the + sign when 
approaching the defect from the right, we see that the bag condition is fulfilled 
when 

= - (43) 
4 + ^72 2 ^^^^ 

thus, g = 2. 

On the other hand, when one approaches the defect from the left, one also 
obtains g = 2 to satisfy the corresponding bag condition. An entirely anal- 
ogous derivation yields the same value for g when the fermions are massive, 
or when more dimensions are considered. For example, in 2 + 1 dimensions, 
and for m = the propagator reduces to the one of [8|. 
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5 Conclusions 



We derived general expressions for the vacuum energy in the presence of 
'defects', whose role is to impose imperfect bag-like boundary conditions on 
parallel planes. Standard bag conditions may be obtained as a particular 
case. For the case of two plates, we obtain more explicit formulae, which 
yield the Casimir energy as an integral over a single variable. That integral 
is both UV and IR finite. 

We show that, at any given distance between the plates, the energy 
reaches a maximum precisely when the bag condition is satisfied. 

Our way of implementing the calculation relies on the calculation of the 
exact fermion propagator in the presence of the defects. This propagator is 
a perfectly well-defined object, in spite of the fact that the eigensystem for 
the associated Dirac Hamiltonian is not well-defined. This ambiguity has the 
consequence of allowing for different parametrizations of the strength of the 
coupling between the Dirac field and the mirrors. In particular, our method 
produces a fermionic propagator which satisfies bag boundary conditions 
when the coupling constant g = 2. That value corresponds, in the setting of 
reference [5] to an infinite coupling constant. 
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